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Abstract 

We report recent progress on understanding the thermalization of the quark-gluon plasma during the early stage in a 
heavy ion collision. The initially high overpopulation in the far-from-equilibrium gluonic matter ("Glasma") is shown 
to play a crucial role. The strongly interacting nature (and thus fast evolution) naturally arises as an emergent property 
of this pre-equilibrium matter where the intrinsic coupling is weak but the highly occupied gluon states coherently 
amplify the scattering. A possible transient Bose-Einstein Condensate is argued to form dynamically on a rather 
general ground. We develop a kinetic approach for describing its evolution toward thermalization, and based on that 
we find approximate scaling solutions as well as numerically study the onset of condensation. 
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1. Introduction 

Thermalization of the quark-gluon plasma is one of the most challenging problems in current heavy ion physics. 
Starting with two colliding nuclei in a form of color glass condensate with high gluon occupation l/a s below saturation 
scale 2s an d following the initial impact, a subsequent strong field evolution stage (likely with instabilities) till about 
the time 1 / Q s is then succeeded by a far-from-equilibrium gluon-dominant matter, the Glasma. The evolution this 
Glasma stage toward a locally equilibrated quark-gluon plasma (QGP) is strongly indicated by phenomenology to be 
reached on the order of a fermi over c time. With these constraints in mind, we explore here a thermalization scenario 
of emergent strongly interacting matter with weak coupling albeit large aggregate of constituents [HE]]. We notice 
recent intensive discussions on this topic from a variety of approaches, see e.g. If3l l4l [5l [6l m [8l [9l ITCH m~l fT^l [T3l IT4l . 

2. Overpopulation as a key feature of the Glasma 

Let us start by considering the somewhat idealized problem of the evolution in a weakly coupled gluon system that 
is initially far from equilibrium and described by the following (Glasma-type) distribution (with coupling a s « 1): 

f(p<QJ = l/a, , f(p>Q*) = Q (1) 

The most salient feature of this initial gluon system, as identified in [ 1 1, is the high overpopulation l/a s , which bears 
a few important consequences by very general arguments and may hold the key of thermalization. 

The initial occupation at a value as high as 1 /a s coherently amplifies scattering and renders the power counting 
in coupling different from normal situation. For example consider the 2 <-> 2 gluon scattering process in the collision 
integral for a transport equation of gluon distribution f(p): while at weak coupling this process contributes at order 
~ d(al) and therefore is "slow" in bringing the system back to equilibrium, in the "Glasma counting" with / ~ l/a s , 
there will be two factors from f(p) and the resulting collision term scales as ~ a\f 2 ~ o(l), despite how small the 
coupling a s may be. Therefore the superficially strongly interacting nature in the Glasma (as implied by fast evolution 
toward equilibrium) could be an emergent property of the weakly coupled albeit highly overpopulated Glasma. 

A highly nontrivial implication of the high overpopulation is that there are so many more gluons in the Glasma 
than in a thermalized plasma for the same amount of energy that a Bose-Einstein Condensation(BEC) has to occur. 
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Let us examine the overpopulation parameter, defined as a dimensionless combination of the particle number density 
and energy density, n e~ 3 ^ 4 . For the Glasma distribution in Eq.Q one has 

n~Ql/a s , e~Qt/a s , (n ~ l/a> /4 (2) 

in sharp contrast to a thermal Bose gas which has ~ 1. More precisely a massless 

Bose gas has (ne~ 3 / 4 ) B = 30 7r7/ i <3 ' > ~ 0.28 << (Bf 4 ' 4 ) G for realistic values a s < 0.3. Therefore the Glasma is 
significantly overpopulated from the outset and all the excessive gluons with the given amount of energy will have 
to be absorbed into a Bose-Einstein Condensate (BEC) if the system evolution is dominated by elastic processes 
(at least over a certain time window). With this analysis on a rather general thermodynamic ground we expect the 
overpopulated Glasma will thermalize to a distribution like f eq (p) - n c (2n) 3 6 3 (p) + l/(e w ?' r - 1) with eventually a 
condensate density parametrically being n c ~ (2 3 /q , s )(1 - \/a l J A ). 

One more important feature of the initial Glasma is that there is only one scale i.e. the saturation scale Q s which 
divides the phase space into two regions, one with / >> 1 and the other with / << 1. One may introduce two 
scales for characterizing a general distribution: a soft scale A s below which the occupation reaches f(p < A s ) ~ 
l/a s >> 1 and a hard cutoff scale A beyond which the occupation is negligible f(p > A) << 1. For initial Glasma 
distribution one has the two scales overlapping A s ~ A ~ <2 S . The thermalization is a process of maximizing the 
entropy (with the given amount of energy). The entropy density for an arbitrary distribution function is given by 
s ~ J^[(l + f)ln{\ + f) — fln(f)]: this implies that with the total energy constrained, it is much more beneficial to 
have as wide as possible a phase space region with / ~ 1 . By this general argument we expect the separation of the 
two scales A s and A in the Glasma along the thermalization process, toward the situation for a very weakly coupled 
thermal gas of gluons with the soft scale Af ~ a s T and the hard scale A th ~ T separated by the coupling a s . 



3. Kinetic equation and scaling solution for Glasma evolution 

With the above general insights about the evolution in overpopulated Glasma, it is tempting to demonstrate these 
more explicitly and quantitatively. To do that we have developed a kinetic approach assuming dominance of 2 <-> 2 
elastic process [12]. In the small-angle approximation one can derive the following transport equation: 



D t f(?) = t(A 2 s A)3 



^f(p) + -[^)f(Mi+f(p)] 



(3) 



P \ A s 

where £ is an order one constant and the two scales A and A s are introduced and defined as: 

We emphasize the full nonlinearity in the /(l + /) terms that arise from the Bosonic nature of gluons and become 
extremely crucial in the highly overpopulated case. One can see that the Glasma distribution implies A, A s ~ <3(1) and 
the collision term C ~ A^A ~ <3(1) in coupling, again in contrast to thermal case with C ~ bio]). 

To describe the thermalization of the Glasma and inspired by the "dropping-out" of coupling, we first discuss 
possible scaling solution for the distribution function f(p) in the static box case. We assume the following scaling 
form characterized by the two scales that evolve in time: 

f(p<A)~^- , f( P >A)~0 (5) 
a s p 

With this distribution the coupling constant entirely drops out from the transport equation Q and the scattering time 
from the collision integral on the RHS scales as f sca ~ A/A;:. To determine the time evolution of A and A s , we need 
two conditions — that the energy must be conserved and that the scattering time shall scale with the time itself, i.e.: 



A A S A 3 
— ~t , e ■ — 

As ffs 

2 



= constant (6) 



The particle number also must be conserved, albeit with a possible component in the condensate: n — n g + n c ~ 
(A s A 2 /a s ) + n c = constant. The condensate plays a vital role with little contribution to energy while unlimited 
capacity to accommodate excessive gluons. Finally with the two conditions we obtain the following scaling solution: 



From this solution, the gluon density n g decreases as ~ (to/t) 1 ^ 7 , and therefore the condensate density is growing with 
time, n c ~ (Ql/a s )[l - (to/f) 1 ' ]. A parametric thermalization time could be identified by the required A s /A ~ a s : 



1 / 1 



7/4 

rth — l-l (8) 

At the same time scale the overpopulation parameter ne~ 3 ^ 4 also reduces to be of order one. 

What would change if one considers the more realistic Glasma with boost-invariant longitudinal expansion? First 
of all the conservation laws will be manifest differently: the total number density will decrease as n ~ noto/t, while 
the time-dependence of energy density depends upon the momentum space anisotropy e ~ %)(to/f) 1+d for a fixed 
anisotropy 5 = Pl/€ (with P L the longitudinal pressure). Along similar line of analysis as before with the new con- 
dition of energy evolution we obtain the following scaling solution in the expanding case: A s ~ Q s (?o/f) (4+(5 ^ 7 , A ~ 
Qs(to/t) (1+26>n - With this solution, we see the gluon number density n g ~ {Ql /a^it^/t)^ 56 ^ 1 , and therefore with 
any 5 > 1/5 the gluon density would drop faster than ~ ?o/f and there will be formation of the condensate, i.e. 
n c ~ (2s/ a s)(fo/0[l _ ( f o/0 <5<5 _1 ^ 7 ]. Similarly a thermalization time scale can be identified through the separation 
i / i \ 7 /( 3 -< 5 > 

of scales to be: ~ -k- \ I . Of course the possibility of maintaining a fixed anisotropy during the Glasma 
evolution is not obvious but quite plausible due to the large scattering rate ~ h\ /A ~ 1/f that is capable of competing 
with the ~ 1/f expansion rate and may reach a dynamical balance: see EQ3] for detailed discussions. 



4. Kinetic evolution toward the onset of condensation 

Of particular significance is to understand how dynamically the condensation occurs starting from an overpopu- 
lated initial condition. With the derived equation such question could be answered by numerically solving it. In 
the case without initial overpopulation, one indeed can show that the system described by this equation evolves all the 
way to a thermal Bose-Einstein distribution which is the proper fixed point of the collision term. In the overpopulated 
case, there is however the complication of the condensate formation. As is well known in atomic BEC literature, one 
has to separately describe the evolution prior to the onset of condensation (with this equation) and the evolution after- 
wards (with a coupled set of two equations explicitly for condensate and regular distribution). Efforts have been made 
in deriving these equations for the coevolution of a condensate and regular distribution [2 1, but here let us focus on the 
pre-BEC stage and investigate how the system approaches the onset of condensation. This stage is solely described 
by the Eq.([3]l and we have numerically solved it for both the static box and the expanding cases. 

With detailed results to be reported in Q, let us highlight the main observation. If the system starts with high 
initial overpopulation, a particle flux in momentum space toward the infrared region will quickly develop and pile up 
particles there. The high occupation number at IR (and thus very fast scattering rate) leads to an almost instantaneous 
local "equilibrium" near the origin p = 0. This local "equilibrium" takes the form: f(p — > 0) — > T*/(p - p*) - 1/2 
with some parameters T* and p* one may tentatively call the local "temperature" and "chemical potential". This form 
is by no means a coincidence — it consists the leading terms in the small p expansion of the Bose distribution which 
is a fixed point of the collision term. With more and more particles being piled up near the origin, the negative local 
"chemical potential" keeps reducing its absolute value and approaches zero i.e. (-p*) — > + . This ultimately marks 
the onset of the condensation. All these have been explicitly seen in the numerical solutions and shown in FigJT] 

We have done extensive numerical studies for varied conditions. For the static box case, the system always reaches 
onset as long as the initial overpopulation parameter ne~ 3 ^ 4 is greater than the thermal value, ne -3 ' 4 > 0.28, despite any 
shape of the initial distribution (e.g. Glasma versus Guassian) or any initial anisotropy. The evolution toward onset 
persists in the expanding case provided enough overpopulation: if starting with isotropic distribution, the critical 
initial overpopulation is shifted mildly to be about ne~ 3 / 4 > 0.40. In the expanding case with initial anisotropy: more 
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Figure 1 : (color online) The time evolution starting with overpopulated initial distribution toward onset of BEC: the left panel shows the distribution 
functions at varied time moments (as blue curves) with red dashed lines showing fitting curves /* = T* /(p - fi*) - 1/2 for the small p region; the 
right panel shows the time dependence of so-extracted "local temperature" T* (upper) and "local chemical potential" fi* (lower). 



initial longitudinal pressure P L > e/3 will shift the critical overpopulation to be smaller, while less initial longitudinal 
pressure Pl < e/3 will shift it to be larger yet only slightly even for large anisotropy. We therefore see that the link 
from initial overpopulation to the onset of condensation, in the present kinetic evolution, is a very robust one. 

Strong evidences for the formation of Bose condensate have been reported for similar thermalization problem 
in the classical-statistical lattice simulation of scalar field theory fl5]|6l|3. The case for non-Abelian gauge theory is 
more complicated and still under investigation fl8] [9] ElQl [HI - in the kinetic approach, the role of inelastic processes 
needs further clarification though there are good arguments for the robustness of condensate even with their presence 
lIU |T5l . We expect many exciting progresses yet to come along this line toward understanding the thermalization. 
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